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We study electron transport through a semiconductor quantum ring with one input and two output 
terminals for an elastic scatterer present within one of the arms of the ring. We demonstrate that 
the scatterer not only introduces asymmetry in the transport probability to the two output leads but 
also reduces the visibility of the Aharonov-Bohm conductance oscillations. This reduction occurs 
in spite of the phase coherence of the elastic scattering and is due to interruption of the electron 
circulation around the ring by the potential defect. The results are in a qualitative agreement with 
a recent experiment by Strambini et al. [Phys. Rev. B 79, 195443 (2009)]. We also indicate that 
the magnetic symmetry of the sum of conductance of both the output leads as obtained in the 
experiment can be understood as resulting from the invariance of backscattering to the input lead 
with respect to the magnetic field orientation. 
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I. INTRODUCTION 

Although studies of electron transport in semiconduc- 
tor rings have a long history!^ the interest in this field 
is sustained by progress of experimental techniques. In 
particular, transport through quantum rings containing 
a number of confined electrons was realized within the 
last decade^ and double concentric quantum rings were 
recently studied^ Moreover, self-interference of electrons 
injected individually into the quantum ring was observed 
with a time-resolved technique^ The Fermi level wave 
functions were probed by conductance measurements for 
the ring potential landscape perturbed by a tip of atomic 
force microscope^ The effect of magnetic forces on quan- 
tum ring conductance was studied experimentally in Ref. 
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In presence of the magnetic forces the electron wave 
function enters both arms of the quantum ring with an 
unequal amplitude £ For two-terminal rings the prefer- 
ential injection of the electron wave function into one 
of the arms of the ring leads to attenuation of the 
Aharonov-Bohm oscillation at high magnetic field3 It 
was demonstrated^ that for rings with three terminals 
at high field in addition to the vanishing oscillation am- 
plitude the magnetic forces produce a distinct imbalance 
of the electron transport probabilities to the two out- 
put leads. Both the high field reduction of the oscilla- 
tion amplitude and the imbalance in the conductance 
of the two output leads were indeed found in the re- 
cent experiment^ However, the experimental data£ differ 
from the theoretical results^ within the range of weak 
magnetic fields, namely: i) the measured conductance of 
one of the output leads significantly exceeds the other 
near B = (Fig. 1 of Ref. Q) and ii) already for low 
magnetic fields the experimental Aharonov-Bohm con- 
ductance oscillations have a low amplitude. The first fea- 
ture suggests that the potential landscape within the ring 
is asymmetric and the second was attributed^ to decohcr- 



ence. The estimated^ coherence length is 320 nm, which 
is surprisingly short - an order of magnitude shorter than 
the estimate for the two-dimensional electron gasii for 
the temperature of 350 mK applied in the experiment^ 
In the present paper we indicate that the observed fea- 
tures of the conductance can also be explained for purely 
coherent transport as resulting from the elastic scattering 
effects which do not randomize the phase but reduce the 
circulation of the electron around the arms of the ring. 
We perform a systematic study of the electron transport 
in a three-terminal ring containing a potential defect. We 
find that only a repulsive and not an attractive scattering 
center may explain the conductance features as seen in 
the experiment. 

The sum of conductance of both the output leads turns 
out^ to be an even function of the magnetic field, which 
is reminiscent of the Onsager symmetry for two-terminal 
devices The sum of the transfer probabilities to the left 
Ti and right T r output leads T = Tj + T r as found in the 
simulation of symmetric ringsi^ is also an even function 
of B, but for evident kinetic reasons which no longer hold 
for asymmetric rings. We demonstrate below that for the 
ring with a defect the kinetics of the electron transfer 
and the electron trajectory are very different for oppo- 
site magnetic field orientations and that the observed^ 
T(B) = T{—B) symmetry is due to the invariance of the 
backscattering with respect to the orientation of the B 
vector. 

Multiterminal rings constitute basic elements for con- 
struction of arrays, which are used in detection of the 
Aharonov-Casher effect^ and are attractive for construc- 
tion of programable quantum gates ^ The three-terminal 
quantum rings were investigated in the context of the 
Kondo density of states^ It was demonstrated that the 
spin-orbit coupling effects in three terminal rings can be 
used for construction of electron spin beam splitters^ 

The current as carried by a steady electron flow at 
the Fermi level can be determined from the Hamiltonian 
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cigcncquation. However, for the purpose of the presen 
study we choose to employ a time dependent approac] 
providing a clear picture of the electron trajectory whicl 
in the present problem appears as quite complex. Wit! 
the wave packet description of the electron motion on 
can approach the time independent monoenergetic limi 
arbitrarily close. 



II. THEORY 

We consider a quantum ring of a radius 155 nm wit! 
three symmetrically attached terminals - see Fig. [IJa) 
The simulations are based on the solution of the time 
dependent Schrodinger equation ih^- = H^> , for th 
Hamiltonian 



600 



H = {-ihV + eA(r)) 2 /2m* 



(1 



in which we apply the Lorcntz gauge A = (—By, 0, 0) an< 
the GaAs effective mass to* = 0.067too. The problem i 
solved with a technique previously used in Refs. ^ n 10 ' 17 
in which the wave function is expanded in a basis 
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of Gaussian functions fj localized around centers (Xj , Yj 



fj(x,y) = Cexp 
ieB , 

+ -yr(x-Xi)(v + Y i ) 



--m^^x-Xtf + iy-Ytf) 



(3 



where C is the normalization constant, ujq determine 
the localization of the basis functions and the imaginar 
term in the exponent introduces the magnetic translatioi 
phase shift that guarantees the gauge invariance, i.e. th 
equivalence of all the applied centers in external magneti 
field. 

The applied choice of centers is shown by the dots in 
Fig. [TJa). The centers are spaced by 22 nm along the 
leads which is close enough to allow for a smooth elec- 
tron flow along the channel provided that the wave vec- 
tor is lower than 0.15 / nm. The electron wave function 
is confined in the direction perpendicular to the axis of 
the leads, and the channel width can be estimated as 
w = 4y/h/m*u)Q. The present modeling of the leads as a 
chain of functions ([3]) limits the simulation to the lowest 
subband. 

The Hamiltonian ([1]) does not explicitly contain any 
confinement potential. In the present model the electron 
confinement results from the localization of the Gaussian 
functions (|3|). Nevertheless, one can try to extract an 
effective confinement potential present within the model 
by considering the Hamiltonian eigenstates obtained in 
basis (|2|). For that purpose we take the eigenequation 
H^ n = and plug expansion Q for the eigen- 

state Subsequently, the weak form of the eigenequa- 
tion is obtained by its projection on the basis elements 
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FIG. 1: (a) The model of the three-terminal quantum ring. 
The dots indicate the positions of the centers of the Gaussian 
basis (0. The green area shows the estimated confinement 
region accessible to the traveling electron (see text) for the 
width of the channel equal to w — 79.2 nm. The effective 
confinement potential is plotted with the blue levels (the low- 
est level corresponds to 10 meV, and the next are spaced by 
40 meV). The colorscale is given in meV. (b) Harmonic os- 
cillator potential (dots) for Tiuq = 2.9 meV and the effective 
potential (solid line) calculated along y — —600 nm line of 
panel (a), (c) The thick curves shows the boundaries of the 
confinement region and the contour plot shows the ground- 
state wave function bound at the junction of the ring to the 
leads. 



([3]), which produces the generalized eigenvalue problem 
Hc„ — E n Sc n , where H and S are the Hamiltonian and 
overlap matrices with analytically intcgrable elements 
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ifi\ H \fj) and s y = (fllfj), respectively The 
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effective potentiai in a point (x, y) is then estimated by 



V(x,y) = 



E n - {-ihV + eA) /2m* y) 



^n{x,y) 



(4) 



The effective potentiai is plotted in Fig. [lja) for Tiwo = 
2.9 meV with the blue contour plot. Additionally a cross 
section of Fig. 0Ja) calculated along the y = —600 nm 
line is plotted in Fig. [IJb) with the solid line. The dots 
in Fig. [TJb) show the harmonic oscillator confinement 
potential m*ui 2 x 2 /2. We can see that the effective con- 
finement potential is consistent with the nominal value 
of cjo applied in the Gaussian basis ©. In Fig. []Ja) 
we marked the channel region in green. The green area 
in the figure was determined as the one in which the 
sum of the Gaussian basis functions exceeds 10% of its 
maximum value, which well agrees with the value of the 
channel width w = 79.2 nm obtained for tkoo = 2.9 meV. 

The junctions of the wire to the ring allow for forma- 
tion of bound electron states. We find three nearly degen- 
erate bound states - a single bound state for a single junc- 
tion - in consistence with the known property of connec- 
tions in the T- wires The binding energy is equal to 0.1 
meV and the wave function of the lowest-energy bound 
state is plotted in Fig. |ljc). For the defect potential 
we use V d (x,y) = W exp [-{{x - X c ) 2 + (y - Y c ) 2 )/R 2 ], 
where (X c , Y c ) are the coordinates of the center of the 
defect, Rd = 28 nm is its radius, and W its height / 
depth. 

As the initial condition for our calculation we take a 
Gaussian wave function entirely localized in the input 
lead (the one below the ring with axis x = 0) localized 
in the direction perpendicular to the axis as the basis 
elements ([3]) but with a larger spread along the lead 

*(r,t = 0) = /,•(*,!/) 

x exp +-m*uj {y - Y 3 f - — — (y - Y 3 ) 2 



x exp(igy), 



(5) 



where the last term in Eq. (JSJ) is introduced to push 
the electron in the direction of the ring with an average 
momentum Tiq. For the applied gauge the kinetic and 
canonical momentum in the y direction are identical and 
the y component of the initial probability density cur- 
rent integrated over the channel is equal to Tiq/m* . The 
Fourier transform of the initial condition along the axis 
of the lead is *(fc) = ^VjllXk exp(-(q - k) 2 /Ak 2 ), and 
AA; is interpreted as the dispersion of the packet in the 
wave vector space. The initial condition |(SJ) is projected 
onto the basis $Z§, and the rest of calculation amounts in 
determining the coefficients c 3 (t) in subsequent moments 
in time. We use the matrix version of the Askar-Cakmak 
scheme^— in form of a system of linear equations for 
c(i + dt), 



Sc(t + dt) = Sc(t - dt) 



lidt 



Hc(t). 



(6) 



We use dt = 0.01 ps. Reduction of the time step below 
this value docs not change the results. 

In the present approach the transfer probabilities Tj 
and T r are determined by the parts of the packet which 
arc transferred to the leads before the end of the simu- 
lation. The simulation is terminated when the electron 
packet completely leaves the ring. We consider the ring as 
empty when it does not contain more than 0.001% of the 
electron charge. Generally, in the time-dependent cal- 
culations the transferred and back-scattered wave pack- 
ets return to the ring after reflection from the ends of 
the channels unless absorbing^ or transparent^ bound- 
ary conditions are used. Application of open bound- 
ary conditions is crucial for approaches using finite dif- 
ference techniques.— £L The present work could be per- 
formed without any open boundary conditions since in 
the present approach one can apply leads of a length 
which is in practice arbitrarily large. A LU decomposi- 
tion of the S matrix for the system of equations ((6|) is per- 
formed before the time stepping. With the decomposed 
overlap matrix the numerical cost of each time step scales 
linearly with the number of centers. For the present cal- 
culation we use in total 5000 Gaussian functions ([3]) with 
the leads as long as 30 fim each. 

In strictly one-dimensional modeling of quantum rings 
with the scattering matrix formalism^ the ring and the 
leads are treated as separate objects with the coupling 
strength described by an appropriate parameter. The 
coupling strength is responsible for the time spent by the 
electron within the ring and the sharpness of the trans- 
fer probability extrcma in function of the wave vector. 
In the present two-dimensional model the ring and the 
leads are modeled as a single object and the ring is es- 
sentially open. Nevertheless the junctions of the ring to 
the leads act like small scattering cavities^ In this pa- 
per the junctions are modeled as right-angle connections. 
Smooth junctions - binding more than a single electron 
state - were recently studied in Ref. [2l|. The type of 
the junction affects the wave vector dependent transfer 
probabilities, but their magnetic-field behavior remains 
qualitatively unchanged. 



III. RESULTS AND DISCUSSION 

We first briefly present the results obtained for an ide- 
ally symmetric configuration (section III. A) to set the 
reference point for the discussion of the transport for a 
defect present in one of the arms of the ring (section 
III.B). 



A. Clean ring 

In order to set the wave vector dispersion parameter 
Afc close enough to the monocncrgctic limit we stud- 
ied the transmission probability through an ideally sym- 
metric ring for B = - see Fig. [5J In the small 
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Afc limit we notice that the variation of T(q) becomes 
more pronounced and the dependence distinctly satu- 
rates. The main features of the saturated T(q) depen- 
dence are well resolved already for the Afc = 1.1 x 1CP 3 
/nm, the value which is used in the rest of the paper. 
Below - unless explicitly stated otherwise - we also as- 
sume q = 0.037 /nm as the average wave vector - the 
value which corresponds to a maximum of T(q) and 
w = 79.2 nm. The kinetic energy of progressive mo- 
tion equals H 2 q 2 /2m* = 0.77 meV, and the dispersion 
of the wave vector assumed covers the energy window 
(h 2 {q-Ak) 2 /2m*,h 2 (q+Ak) 2 /2m*) = (0.73,0.82) meV. 
The second subband is by tiujo = 2.9 meV above the low- 
est one - far above the electron kinetic energy - which 
is consistent with the neglect of the scattering to higher 
subbands assumed in the present approach. 
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q[ 1 /nm ] 
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q[ 1 / nm ] 

FIG. 2: Electron transfer probability to the left T; and right 
T r leads for B = in function of the average wave vector q 
for a number of Afc values for w = 79.2 nm. 

The effect of the magnetic field on the wave- vector re- 
solved transfer probability through the symmetric ring is 
illustrated in Fig. [3J We notice that the external mag- 
netic field introduces asymmetry of the transfer to the 
left and right output leads. For B > 0, TJ increases at 
the expense of the T r , which is a direct consequence of 
the magnetic forces which preferentially inject the elec- 
tron to the left arm of the ring and then eject it to the 
left output lead. For higher magnetic fields T; and T r 
become less strongly dependent on q. 

The dependence of the transfer probability on B is 
shown in Fig. HJa). The transmission probabilities 2], T r 
as well as their sum T undergo oscillations that are due to 
the Aharonov-Bohm effect. The period of the oscillations 
of T is equal to 0.055 T, in agreement with the nominal 
value of the magnetic field, which corresponds to the flux 
quantum $o = e /h for the ring of radius 155 nm. The 
distinct decrease of the oscillations amplitude for larger 
B is another consequence of the magnetic forces. For 
higher fields most of the electron packet is injected into 



FIG. 3: Electron transfer probability to the left Tj and right 
T r lead for B = 0, 0.3 and 0.6 T in function of the average 
wave vector q for Afc = 1.1 x 10 -3 /nm. 

one of the arms of the ring, which subsequently intro- 
duces an imbalance in the parts of the packet that meet 
and interfere near the output leads. In consequence the 
Aharonov-Bohm interference is less pronounced. 

In order to demonstrate the effect of the channel width 
we presented in Fig. 0Jb) and (c) the transfer probabili- 
ties for wider channels, namely for w = 90.5 nm (b) and 
w = 113.13 nm (c) (in Fig. [4] and everywhere else in 
this paper we assume w = 79.2 nm). The strength of the 
magnetic deflection of the electron trajectories depends 
on the ratio of the Larmor radius to the channel width 
- for wider channels there is more space for the mag- 
netic deflection, and consequently the Aharonov-Bohm 
oscillations vanish faster for larger w. 

In Fig. [3] we also notice, that even for the symmetric 
ring the extrema of Ti and T r are shifted off B = (in 

other words -gg^- \ b=o 0) onr y the sum of the trans- 
mission probabilities T(B) = Ti(B) + T r (B) is an even 
function of B. For the symmetric ring one has 

T l (B)=T r (-B), (7) 

which suffices to explain the T(B) = T(—B) symmetry 
observed in Fig. [3J However, for asymmetric rings rela- 
tion j7]) no longer holds, although the transfer probability 
remains symmetric in B - see below. 

B. Ring with a defect 

Before inserting the Gaussian defect to the ring we 
first study its scattering properties for a straight channel. 
Fig. [5] shows the transmission probability as function of 
the height of the defect W for the average wave vector 
q = 0.037 / nm. For the attractive defect (W < 0) the 
transmission probability is close to 1, and the defect acts 
like a phase shifter. For W > the defect is a more 
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FIG. 4: Electron transfer probability to the left Ti and right 
T r leads as well as their sum (shifted up for clarity by 0.5) for 
pure ring and channels of width w = 79.2 nm (a) w = 90.5 
nm (b) and w = 113.13 nm (c). 



effective scatterer with transmission probability as low 
as 0.08% for W = 5 mcV. 

Next, we place a defect of height W = 3 meV in the 
left arm of the ring between the input and the left output 
lead - see the red circle in Fig. [6j which also shows the 
amplitude of the electron wave function for B = ±0.6 T 




-5 5 
W [meV] 



FIG. 5: Electron transfer probability for a straight chan- 
nel containing a Gaussian defect of radius Rd = 28 nm and 
height/depth W (negative values of W correspond to poten- 
tial well), for B = 0. 



for several moments in time. For the negative B most 
of the electron wave function is injected to the right arm 
and next to the right output lead, as for a pure quantum 
ring. The dominant trajectory for this field is drawn 
schematically in Fig. [7Ja). More complex is the transport 
for the positive magnetic field [the lower panel of Fig. [6] 
and Fig. E£b)]: the electron is first injected into the 
left arm, then it is nearly completely reflected by the 
defect. The electron velocity is inverted and the Lorentz 
force - still tending to deflect the trajectory to the left 
- keeps the electron within the ring as it passes near 
the input lead and then both output leads. The electron 
is subsequently reflected again from the defect, this time 
from its other side. Only after the second scattering event 
the magnetic force pushes the electron to the left output 
ring. In consequence most of the probability density goes 
to the left output lead just like for the pure ring. In 
consistence with the schematic drawing of the dominant 
trajectory of Fig. [T](b) in Fig. [5] for B = 0.6 T we 
first observe an increased probability amplitude below 
the defect (48.6, 54, 64.8 ps) and then above it (75.7 and 
86.5 ps). 

Fig. [8] shows the parts of the electron packet within 
the ring and in the leads calculated for the simulation 
presented in Fig. [Sj For the positive magnetic field the 
electron trajectory has a larger length [Fig. [7jb)] and 
the electron packet stays longer within the ring. In the 
large t limit we see in Fig. [5] that relation ([7]) no longer 
holds. However, the part of the electron packet in the 
input lead (the incoming and backscattered parts of the 
packet), is for any t exactly the same for both magnetic 
field orientations. The independence of the backscatter- 
ing of the magnetic field orientation can be understood 
as due to the fact that the backscattered trajectories are 
identical for ±B ^ As a result of this invariance, the re- 
lation T(B) = T(—B) holds in spite of the very different 
kinetics of the electron transfer through the asymmetric 
ring for B = ±0.6 T. 

The oscillations of the transfer probabilities in func- 



6 




x [nm] 



FIG. 6: Amplitude of the wave function calculated for B = —0.6 T (the upper row of plots) and B = 0.6 T (the lower row) at 
subsequent moments in time (given on top of the figure). The circle is centered at the position of the repulsive defect (W = 3 
meV). 
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FIG. 7: Schematic drawing of the dominant trajectory in the 
electron transport presented in Fig. [S]for B = —0.6 T (a) and 
B — 0.6 T (b). 
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FIG. 8: Part of the wave packet inside (red curves) the ring, 
in the input lead (black curves), in the left (blue curves) and 
right (green curves) output leads. Results for B — 0.6 T are 
plotted with dashed curves and for B = —0.6 T with solid 
ones. The data correspond to the wave function snapshots 
presented in Fig. [5] 

tion of B arc plotted in Fig. [9] for attractive (a-d) and 
repulsive (e-h) defects of different height / depth. The 
results can be compared with the ones given for the pure 
ring in Fig. Hla), and the parameters considered in Figs. 



HI [7] and [5] are applied in Fig. dig). The presence of 
the attractive defect, which (only) shifts the phase of the 
part of the wave function passing through the left arm, 
changes the local maximum of T for B = [Fig. HJa) and 
HJd)] into a local minimum [Fig. Hfb) and (c)]. However, 
the average values of 7] and T r within the range of small 
magnetic fields B G [—0.2 T, 0.2 T] remain very simi- 
lar. On the contrary, the repulsive defect [Fig. SJc-h)] 
leads to a pronounced difference in T; and T r values for 
B ~ 0. Moreover, since the circulation of the electron 
around the ring is stopped by the repulsive defect (see 
Fig. [5]), its presence drastically reduces the amplitude of 
the Aharonov-Bohm oscillations. 

The results of Fig. [UJb) are in a good agreement with 
the experimental results for conductance given in Fig. 1 
of Ref. [H within the range of the magnetic fields pre- 
sented therein: i) the amplitude of the oscillation is small, 
ii) for low magnetic fields T r distinctly exceeds T; , iii) 7] 
dominates for B > and T r for B < iv) the oscillations 
of the transfer probabilities vanish at higher magnetic 
field, v) the overall transfer probability T stays symmet- 
ric in B, vi) the envelope of the T oscillations possesses 
a pronounced minimum near B = 0. 

In order to establish which of the above features are 
independent of the wave vector we plotted in Fig. [TOlthe 
transfer probabilities for B = in function of q for fixed 
value of the height of the defect W = 3 meV, and in Fig. 
1111 the transfer probabilities in function of the magnetic 
field for several fixed values of q. In Fig. [10] we see that 
the defect present in the left arm of the ring does not 
necessarily imply Ti <T r , for B = 0. Results of Fig. [Tj] 
(a) and (d) were calculated for local maxima of T/. The 
crossings between XJ and T r appears here for the negative 
magnetic field instead of the positive B as in Fig. [5fe- 
g) . The overall transfer probability has a maximum near 
B = for q = 0.032 / nm, and for q = 0.048 / nm the 
value of T oscillates around an average constant value. 
For q = 0.042/nm [Fig. FITTc)] near the minimum of 
T(q) (Fig. [TO]) Ti and T r have similar values on a longer 
range of B between and 0.2 T. The presented results 
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FIG. 10: Electron transler probability to the left Ti and right 
T r leads for the ring with a repulsive defect of height 3 meV 
and B = in function of the wave vector q. 



indicate that of the above list of features only i) , iii) and 
v) are characteristic to the ring with a strong scatter- 
ing repulsive defect, and the rest is g-dependent. The 
Fermi wave vector in the experimental samples is deter- 
mined by the density of the two-dimensional electron gas 
in the electron reservoirs. In principle the electron den- 
sity is fixed at the sample formation stage by the dopants 
concentration within the AlGaAs barrier, but the wave 
vector should be at least to an extent tunable by the 
voltages applied to the electrodes in a gated sample. 

In order to comment on the location of the defect that 
we consider here, we indicate that a defect in between 



the input and the right output lead produces the tran- 
sition spectra with inverted value of the magnetic field 
(B — > —B) and a defect placed exactly in the center of 
the arm between the output leads produces symmetric 
spectra only with oscillations amplitude that is reduced 
for W > 0. For B = the exact position of the defect 
within the same section of the ring influences the posi- 
tion of the transfer probability peaks on the wave vector 
scale, but otherwise no qualitative difference is found in 
the transfer characteristics as function of the external 
magnetic field. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the electron transport through a 
three-terminal quantum ring containing clastic scatterers 
using a time-dependent approach. The presented study 
indicates that elastic scattering may be a significant rea- 
son of the low amplitude of Aharonov-Bohm oscillations 
simultaneously explaining the low-field asymmetry of the 
conductance to both the output leads as observed in a 
recent experiment. Low visibility of the oscillations in- 
troduced by the elastic scattering is not due to the phase 
randomization but to a hindered circulation of the elec- 
tron around the ring. Both decohcrence and intersub- 
band scattering^ that were neglected in the presented 
study should also reduce the visibility of the Aharonov- 
Bohm oscillations but by themselves they cannot intro- 
duce the strong asymmetry effects that are distinct in 
the experimental results. We also found that the poten- 
tial defect present within the ring affects the properties 
of the transmission probability only in the low magnetic 
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FIG. 11: Electron transfer probability for the repulsive defect W — 3 meV for average wave vector values q — 0.032/nm (a), 
q = 0.038/nm (b), q = 0.042/nm (c) and q = 0.048/nm (d) in function of the magnetic field. 



field. The high magnetic field limit is left unchanged: the 
Aharonov-Bohm oscillations are reduced by the electron 
injection imbalance due to the magnetic forces and the 
conductance of one of the leads increases at the expense 
of the other as in the case of a clean quantum ring. 
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